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Abstract 

Low's well known action principle for the Maxwell- Vlasov equations of 
ideal plasma dynamics was originally expressed in terms of a mixture of Eule- 
rian and Lagrangian variables. By imposing suitable constraints on the vari- 
ations and analyzing invariance properties of the Lagrangian, as one does for 
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the Euler equations for the rigid body and ideal fluids, we first transform this 
action principle into purely Eulerian variables. Hamilton's principle for the 
Eulerian description of Low's action principle then casts the Maxwell- Vlasov 
equations into Euler-Poincare form for right invariant motion on the diffeo- 
morphism group of position- velocity phase space, R 6 . Legendre transforming 
the Eulerian form of Low's action principle produces the Hamiltonian formula- 
tion of these equations in the Eulerian description. Since it arises from Euler- 
Poincare equations, this Hamiltonian formulation can be written in terms of 
a Poisson structure that contains the Lie-Poisson bracket on the dual of a 
semidirect product Lie algebra. Because of degeneracies in the Lagrangian, 
the Legendre transform is dealt with using the Dirac theory of constraints. 
Another Maxwell- Vlasov Poisson structure is known, whose ingredients are 
the Lie-Poisson bracket on the dual of the Lie algebra of symplectomorphisms 
of phase space and the Born-Infeld brackets for the Maxwell field. We dis- 
cuss the relationship between these two Hamiltonian formulations. We also 
discuss the general Kelvin-Noether theorem for Euler-Poincare equations and 
its meaning in the plasma context. 
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1 Introduction 

Reduction of action principles. Due to their wide applicability, the Maxwell- 
Vlasov equations of ideal plasma dynamics have been studied extensively. Low 
[1958] wrote down an action principle for them in preparation for studying stability 
of plasma equilibria. Low's action principle is expressed in terms of a mixture of 
Lagrangian particle variables and Eulerian field variables. 
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Following the initiative of Arnold [1966] and its later developments (see Marsdcn 
and Ratiu [1994] for background), we start with a purely Lagrangian description of 
the plasma and investigate the invariance properties of the corresponding action. 
Using this set up and recent developments in the theory of the Euler-Poincare 
equations (Poincare [1901b]) due to Holm, Marsden and Ratiu [1997], we are able 
to cast Low's action principle into a purely Eulerian description. 

In this paper, we start with the standard form of Hamilton's variational principle 
(in the Lagrangian representation) and derive the new Eulerian action principle by 
a systematic reduction process, much as one does in the corresponding derivation of 
Poisson brackets in the Hamiltonian formulation of the Maxwell- Vlasov equations 
starting with the standard canonical brackets and proceeding by symmetry reduction 
(as in Marsdcn and Weinstein [1982]). In particular, the Eulerian action principle 
we obtain in this way is different from the ones found in Ye and Morrison [1992] by 
ad hoc procedures. We also mention that the method of reduction of variational 
principles we develop naturally justifies constraints on the variations of the so called 
"Lin constraint" form, well known in fluid mechanics. 

The methods of this paper are based on reduction of variational principles, 
that is, on Lagrangian reduction (see Cendra et al. [1986, 1987] and Marsden and 
Scheurle [1993a,b]). These methods have also been useful for systems with non- 
holonomic constraints. This has been demonstrated in the work of Bloch, Krish- 
naprasad, Marsden and Murray [1996], who derived the reduced Lagrange d'Alem- 
bert equations for nonholonomic systems, which also have a constrained variational 
structure. The methods of the present paper should enhance the applicability of 
the Lagrangian reduction techniques for even wider classes of continuum systems. 

Passage to the Hamiltonian formulation. The Hamiltonian structure and 
nonlinear stability properties of the equilibrium solutions for the Maxwell- Vlasov 
system have been thoroughly explored. Some of the key references are Iwinski and 
Turski [1976], Morrison [1980], Marsden and Weinstein [1982] and Holm, Marsden, 
Ratiu and Weinstein [1985]. See also the introduction and bibliography of Marsden, 
Weinstein et al. [1983] for a guide to the history and literature of this subject. 

In our approach, Lagrangian reduction leads to the Euler-Poincare form of the 
equations, which is still in the Lagrangian formulation. Using this set up, one may 
pass from the Lagrangian to the Hamiltonian formulation of the Maxwell- Vlasov 
equations by Legcndrc transforming the action principle in the Eulerian description 
at either the level of the group variables (the level that keeps track of the particle 
positions), or at the level of the Lie algebra variables. One must be cautious in this 
procedure because the relevant Hamiltonian and Lagrangian are degenerate. We 
deal with this degeneracy by using a version of the Dirac theory of constraints. 

Legendre transforming at the group level leads to a canonical Hamiltonian for- 
mulation and the latter leads to a new Hamiltonian formulation of the Maxwell- 
Vlasov equations in terms of a Poisson structure containing the Lie-Poisson bracket 
on the dual of a semidirect product Lie algebra. This new formulation leads us 
naturally to the starting point for Hamiltonian reduction used by Marsden and 
Weinstein [1982] (see also Morrison [1980] and Kaufman and Dewar [1984]). 
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Stability and asymptotics. The new Hamiltonian formulation of the Maxwell- 
Vlasov system places these equations into a framework in which one can use the 
energy-momentum and energy-Casimir methods for studying nonlinear stability 
properties of their relative equilibrium solutions. This is directly in line with Low's 
intended program, since the study of stability was Low's original motivation for 
writing his action principle. Sample references in this direction are Holm, Mars- 
den, Weinstein and Ratiu [1985], Morrison [1987], Morrison and Pfirsch [1990], Wan 
[1990], Batt and Rein [1993] and Batt, Morrison and Rein [1995]. Other historical 
references for the Lagrangian approach to the Maxwell- Vlasov equations include 
Sturrock [1958], Galloway and Kim [1971] and Dewar [1972]. 

The Eulcrian formulation of Low's action principle also casts it into a form that 
is amenable to asymptotic expansions and creation of approximate theories (such 
as guiding center theories) possessing the same mathematical structure arising from 
the Euler-Poincare setting. See, for example, Holm [1996] for applications of this 
approach of Hamilton's principle asymptotics in geophysical fluid dynamics. 

Comments on the Maxwell- Vlasov system. The rest of this paper will be 
concerned with variational principles for the Maxwell- Vlasov system of equations for 
the dynamics of an ideal plasma. These equations have a long history dating back at 
least to Jeans [1902], who used them in a simpler form known as the Poisson- Vlasov 
system to study structure formation on stellar and galactic scales. Even before 
Jeans, Poincare [1890, 1901a] had investigated the stability of equilibrium solutions 
of the Poisson- Vlasov system for the purpose of determining the stability conditions 
for stcllcr configurations. The history of the efforts to establish stellar stability 
conditions using the Poisson- Vlasov system is summarized by Chandrasekhar [1977] . 
The Poisson- Vlasov system is also used to describe the self-consistent dynamics of 
an electrostatic collisionless plasma, whereas the Maxwell- Vlasov system is used 
to describe the dynamics of a collisionless plasma evolving self-consistently in an 
electromagnetic field. 

Organization of the paper. The paper is organized as follows. Section || in- 
troduces the Maxwell- Vlasov equations. In section || we state the Euler-Poincare 
theorem for Lagrangians depending on parameters along with the associated Kelvin- 
Noether theorem. This general theorem plays a key role in our analysis. Section 
H reformulates these equations in a purely Eulerian form and shows how they sat- 
isfy the Euler-Poincare theorem. The following section reviews some aspects of the 
Legendre transformation for degenerate Lagrangians. Section |J reprises Low's ac- 
tion principle for the Maxwell- Vlasov equations. Section ^ casts the Euler-Poincare 
formulation of the Maxwell- Vlasov equations into Hamiltonian form possessing a 
Poisson structure that contains a Lie- Poisson bracket. In Section ||we summarize 
our conclusions. 
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2 The Maxwell- Vlasov equations 

The Maxwell- Vlasov system of equations describes the single particle distribution 
for a set of charged particles of one species moving self-consistently in an elec- 
tromagnetic field. In this description, the Boltzmann function /(x, v, t) is viewed 
as the instantaneous probability density function for the particle distribution, i.e., 
given a region £1 of phase space, the probability of finding a particle in that region 
is 



/ dxdvf(x,v,t), 
Jn 



(2.1) 



where x and v are the current positions and velocities of the plasma particles. Thus, 
if the phase-space domain is the whole (x,v) space, the value of this integral at 
a certain time t is normalized to unity. 

As is customary, we assume that the particles of the plasma obey dynamical 
equations and that the plasma density / is advected as a scalar along the particle 
trajectories in phase space, i.e., 

dj_ 
dt 

In this equation, and in the sequel, an overdot refers to a time derivative along a 
phase space trajectory, and V x and V v denote the gradient operators with respect 
to position and velocity respectively. For pressureless motion in the electromagnetic 
field of the charged particle distribution, the acceleration of a particle is given by 



x-V x / + v- V v / = 0. 



(2.2) 



x = 



dA 

V x $ + — - v x (V x x A) 



(2.3) 



where (q/m) denotes the charge to mass ratio of an individual particle, $ is the elec- 
tric potential, and A is the magnetic vector potential. Substituting this expression 
for v in equation (2.2) yields 



-rrr + v • V x / 

at m 



dA 

V x $ + — - v x (V x x A) 
at 



V v / = 0. 



(2.4) 



This is the Vlasov equation (also called the collisionlcss Boltzmann, or Jeans 
equation). The system is completed by the Maxwell equations with sources: 



V x • E = p, 



T3 dE 



(2.5) 



where E and B are the electric and magnetic field variables respectively, p is the 
charge density and j is the current density. These quantities are expressed in terms 
of the Boltzmann function / and the Maxwell scalar and vector potentials $ and 
A by: 

dA 

~&t' 



E 



-v x $ 



B = V x x A, 



p(x,t) = q 



J dv /(x,v,t), j(x,t)=q J (fvv/(x,v,t). 



(2.6) 
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By their definitions, E and B satisfy the kinematic Maxwell equations 

V x -B = 0, V x xE = -— . (2.7) 



Equations (2.4) - (2/7) comprise the Maxwell- Vlasov equations. When A is 
absent, the field is electrostatic and one obtains the Poisson- Vlasov equations. The 
Poisson- Vlasov system can also be used to describe a self gravitating collisionless 
fluid, and so it forms a model for the evolution of galactic dynamics, see, e.g., Binney 
and Tremaine [1987]. 



Note that the integral in (2.1) is independent of time (as the region and the func- 
tion / evolve), since the vector field defining the motion of particles (see equation 
( [2.3| )) is divergence free with respect to the standard volume element on velocity 
phase space. Thus, one may interpret / either as a density or as a scalar. For our 
purposes later, we will need to be careful with the distinction, since the volume 
preserving nature of the flow of particles will be a consequence of our variational 
principle and will not be imposed at the outset. 



3 The Euler-Poincare equations, Semidirect Prod- 
ucts, and Kelvin's Theorem 

The general Euler-Poincare equations. Here we recall from Holm, Marsden 
and Ratiu [1997] the general form of the Euler-Poincare equations and their asso- 
ciated Kelvin-Noether theorem. In the next section, we will immediately specialize 
these statements for a general invariance group G to the case of plasmas when 
G is the diffeomorphism group, Diff(TR 3 ). We shall state the general theorem for 
right actions and right invariant Lagrangians, which is appropriate for the Maxwell- 
Vlasov situation. The notation is as follows. 

• There is a right representation of the Lie group G on the vector space V and 
G acts in the natural way from the right on TG x V*: (v g , a)h = (v g h, ah). 

• p v : g — > V is the linear map given by the corresponding right action of the 
Lie algebra on V: p v {0 — v£, and p* : V* — > g* is its dual. The g-action on 
g* and V* is defined to be minus the dual map of the g-action on g and V 
respectively and is denoted by ^£ and a£ for £ G g, p 6 g*, and a G V* . For 
v G V and a G V* , it will be convenient to write: 

voa = pla i.e., (v o a, t) = (a, u£) = - («, a£) , 

for all £ G g. Note that v o a G g*. 

• Let Q be a manifold on which G acts trivally and assume that we have a 
function L : TG x TQ xf-^l which is right G-invariant. 

• In particular, if ao G V* , define the Lagrangian L ao : TG x TQ — > K by 

) = L(vg,u q ,a ). Then L ao is right invariant under the lift to TG x 
TQ of the right action of G ao on G x Q. 
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• Right G-invariance of L permits us to define I : g x T Q x V* ^ R by 

Itygg^iUq.ag' 1 ) = L(v g ,u q ,a). 

Conversely, this relation defines for any I : g x TQ x V* R a right G— 
invariant function L : TG x TQ x V* -> R. 

• For a curve g(t) € G, let := <7(<)g(i) _1 and define the curve a(t) as 
the unique solution of the linear differential equation with time dependent 
coefficients a(t) — —a(t)£(t) with initial condition a(0) = do- The solution 
can be equivalently written as a(t) = ao g{t) . 



Theorem 3.1 The following are equivalent: 
i Hamilton's variational principle holds: 

L ao (g(t),g(t),q{t),q{t))dt = , (3.1) 

for variations of g and q with fixed endpoints. 
ii (g(t),q(t)) satisfies the Euler-Lagrange equations for L ag on G x Q. 
iii The constrained variational principle^ 

l(Z(t),q(t),q(t),a(t))dt = 0, (3.2) 

holds on q x Q, upon using variations of the form 

5£=^- ad^n = ^ - £,77], 5a = -an, (3.3) 

where n(t) 6 g vanishes at the endpoints and Sq(t) is unrestricted except for 
vanishing at the endpoints. 

iv The following system of Euler-P oincare equations (with a parameter) coupled 
with Euler- Lagrange equations holds on g x TQ x V* : 

8 61 , m 5l SI 

mrr-^si + ra^' (3 - 4) 

and 

d dl dl 



1 Strictly speaking this is not a variational principle because of the constraints imposed on 
the variations. Rather, this principle is more like the Lagrange d'Alembert principle used in 
nonholonomic mechanics. 
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The strategy of the proof is simple: one just determines the form of the variations 
on the reduced space g x Q x V* that are induced by variations on the unreduced 
space TG x TQ and includes the relation of a(t) to a . One then carries the 
variational principle to the quotient. See Holm, Marsden and Ratiu [1997] for 
details. Here we have included the extra factor of Q which is needed in the present 
application; this will be the space of potentials for the Maxwell field. This extra 
factor does not substantively alter the arguments. 



The Kelvin-Noether Theorem. We start with a Lagrangian L ao depending 
on a parameter clq £ V* as above and introduce a manifold C on which G acts. 
We assume this is also a right action and suppose we have an equivariant map 
IC:CxV* ^ 0**. 

In the case of continuum theories, the space C is chosen to be a loop space and 
(/C(c, a), [/,) for c £ C and fj, £ g* will be a circulation. This class of examples also 
shows why we do not want to identify the double dual g** with q. 

Define the Kelvin-Noether quantity / : C x g x TQ x -> R by 

I(c,£,q,q,a) = (/C(c,a), ^(€,q,q,a)) . (3.6) 



Theorem 3.2 (Kelvin-Noether) Fixing cq £ C, let , q(t) , q(t) , a(t) satisfy the 
Euler-P oincare equations and define g(t) to be the solution of g(t) = £(t)g(t) and, 
say,g(0) = e. Let c(t) = ff(t)" 1 cq and I(t) = I(c(t),£(t),q(t),q(t),a(t)). Then 

j t I(t) = ^fC{c(t),a(t)),^o a y (3.7) 

The proof of this theorem is relatively straightforward; we refer to Holm, Mars- 
den and Ratiu [1997]. We shall express the relation ([Tt]) explicitly for Maxwell- 
Vlasov plasmas at the end of section R. 



4 An action for the Maxwell- Vlasov equations 

A typical element of TR 3 =M 3 xl 3 will be denoted z = (x, v). We let n s : TR 3 -> 
R 3 and ir v : TR 3 — > R 3 be the projections ir s (z) = x and ir v (z) = v onto the first 
and second factors, respectively. 



Spaces of fields. We let Diff(TR 3 ) denote the group of C°°- diffeomophisms 
from TR 3 onto itself. An element ip £ Diff(TR 3 ) maps plasma particles having 
initial position and velocity (xo,Vo) to their current position and velocity (x, v) = 
y>(xo,Vo). This is the particle evolution map. We shall sometimes abbreviate 
(xo,vo) = zo, (x, v) = z, etc. The spatial components of ^(xo,vo) are written as 
x(xo, vo) and the velocity components as v(xo, vo). We shall also use the following 
notation: 

• V = C°°(R 3 ,R) is the space of electric potentials $(x); 
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• A is the space of magnetic potentials A(x); 

• T = C°°(T]R 3 ,R) is the space of plasma densities /(x, v); 

• To = C^°(TR 3 , M) is the space of plasma densities with compact support; 

• T> = Cq°(R 3 ,R) is a space of test functions, denoted <p(x). 

The test functions <p(x) are used to localize the variational principle. Thus, once 
one obtains Euler-Lagrange equations depending on fo and ipo, if their validity can 
be naturally extended for any fo and ipo, which will happen in our case, then we 
shall consider those extended equations to be the Euler-Lagrange equations of the 
system. We will usually be interested in the Euler-Lagrange equations for fo > 
and ipo = 1. 

The Lagrangian and the action. For each choice of the initial plasma distri- 
bution function fo and the test function ipo, we define the Lagrangian 

Lf 0lV0 {il>,ip,®,®,A.,A.) = J rfx dv /o(x , v ) Qm|x(x , v )| 2 

+ im|x(x ,vo) - v(x ,v )| 2 (4.1) 

+ <?x(x , v ) • A(x(x , v )) - q$(x(x , v )) 

+ \ J dr^r) (|V r $ + ^(r)| 2 - |V r x A(r)| 2 ) . 

This Lagrangian is the natural generalization of that for an iV-particle system, 
with terms corresponding to kinetic energy, electric and magnetic field energy, the 
usual magnetic coupling term with coupling constant q (the electric charge), and a 
constraint that ties the Eulerian fluid velocity v to x, the material derivative of the 
Lagrangian particle trajectory. Here x and v are Lagrangian phase space variables, 
while A and <& are Eulerian field variables. Thus, there should be no confusion 
created by the slight abuse of notation in abbreviating dA/dt and d^/dt as $ and 
A, respectively, in the arguments of the Lagrangian. This Lagrangian is inspired 
by Low [1958]. However, we have added the term 

im|x(x ,v ) - v(x ,v )| 2 , 

which allows v to be varied independently in the variational treatment. 
Consider the action 

6 = J dt Lf 0iVo (i[), tp, <£", A, A) , 

defined on the family of curves {tp(t), $(t), A(i)) satisfying the usual fixed-endpoint 
conditions (ip(ti), $>(U), A(tj)) = (^,$j,Aj), i = 1,2. One now applies the stan- 
dard techniques of the calculus of variations. In particular, integration by parts can 
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be performed since /o and ipo have compact support. Moreover, once the Euler- 
Lagrange equations have been obtained, their validity can be easily extended in a 
natural way for / > and (p = I. 



Derivation of the equations. To write the equations of motion, we need some 
additional notation. Consider the evolution map ^(xojVo) = (x, v) so that ipt 
relates the initial positions and velocities of fluid particles to their positions and 
velocities at time t. Let u be the corresponding vector field: 



u(x,v) := ip t oip t 1 (x,v) =: x 



8 



d 



<9x 9v ' 



so the components of u are (x, v) . Recall that the transport of fo as a scalar is 
given by /(x, v, t) — / o tp^(x, v), which satisfies 



0. 



(4.2) 



where V z = (V X ,V V ) is the six dimensional gradient operator in (x, v) space. 
Let Jif, be the Jacobian determinant of the mapping ip 6 Diff(TM 3 ), that is, the 
determinant of the Jacobian matrix <9(x, v)/<9(xo, vq). 
Define F(x, v,i) to be /o, transported as a density: 



F(x(x , v ), v(x , v ),t)Jv( x o, v ) = /o(xo,v ), 



so that 



dF 
~dt 



V z • (Fu) = 0. 



(4.3) 



Taking variations in our Lagrangian (4.1) and making use of the preceding equation 
for F, we obtain the following equations (taking ipo = 1) 



<5x 
<5v 

SA 



mx + m(x — v) 



x- v = 0, 
dA 



OA 

q— + qxx (V x x A), 



dt 



dt 



= -q 



d-vF(x, v, t), 

d ( dA\ f 

V x x (V x x A) = -— I V x $ + -^rj +q dvvF(x,v,i). 



(4.4) 



The second equation in (4.4) treats the Eulerian velocity v as a Lagrange multi- 
plier, and ties its value to the fluid velocity x, hence v = x as well. The first 
two variational equations in the set (14) provide the desired relation for parti- 
cle acceleration and the last two equations are the Maxwell equations with source 
terms. Thus, Hamilton's principle with Low's action provides the equations for 
self-consistent particle motion in an electromagnetic field, as required, and the de- 
scription is completed by substituting 



m 



V x $ 



dA 

~dt 



v x (V x x A) 
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for the components of u in the transport equation (4.2) to give the Vlasov equation 



5 The Maxwell-Vlasov system as Euler-Poincare 
equations 

We will now specialize the general Euler-Poincare the orem to the case of plasmas. 
The Lagrangian Lf 0tVo (tp, ip, $, $, A, A) in equation (11) has a right Diff(TIR 3 )- 
symmetry. Let r\ € Diff(TR 3 ), F £ T and define the action of r\ on F by Frj — 
(F o i])J v where, as above, J v is the Jacobian determinant of rj. 
The symmetry of L/ Oj¥ , o (-0, ip, $, A, A) is the property 

Lfov,vo (V"7j V"7; ^> A, A) = Lf o ipo (ip,ip,$,$, A, A), 
for all i] e Diff(TlR 3 ). 



Ingredients for Euler-Poincare. Now we apply the general Euler-Poincare the- 
orem 3.1, taking G = Diff(T]R 3 ) and Q = V x A and the parameter eto = /o- As 



we have explained before, tpo is an auxiliary quantity that will ultimately take the 
value unity. In the general Euler-Poincare theorem 3.1 we take 



ou = — - ad u w, 
at 



Sa = — £ v 



(5.1) 



where w G g is a vector field on TR 3 , £ w is the Lie derivative and ad u w = — [u, w] 
defines ad u w i n terms of the Lie bracket of vector fields, [u, w]. The Euler-Poincare 
equations (3.4) are 



d SI w 61 SI 

at Su Su Sa 



(5.2) 



where ad* is the dual of ad u and Sl/5u is a 1-form density. The 1-form density 
{81 /Sa) o a is defined by 



/ 51 



= I — ■ -t w a. 
Sa 



(5.3) 



When the quantities a are tensor fields, 81 /Sa will be elements of the dual space 
under the natural pairing. 



We shall apply this result to obtain the Maxwell-Vlasov system (2.4)-(2.7) as 



Euler-Poincare equations. We begin by recording a formula that will be needed later. 
Let u, w be two elements of g, the Lie algebra of vector fields for the diffeomorphism 
group on a manifold M. Choose the 1-form density c 6 g* , and let the pairing 
(c,u}:g*xg^M be given by: 



c, u) = / dzc ■ u = / i 
Jm J m 



dz Cj 



(5.4) 
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where Cj and u J , j — 1 . . . n, are components of c and u in R™ and dz is the volume 
form on M. Then we can write the desired formula, 



(ad* c, w) = J dz ad* c • w = / dzc ■ ad u w 



' /zr - " 7777 ""'7777 



= j dzw 1 (cj^r + CiiV -u) + (u- v)c, 
= (£ u c,mv). 



(5.5) 



Here ,£ u c is the Lie derivative of the the 1-form density c with respect to the vector 
field u, z J is the coordinate chart and Cj ,u J 1 w J are the components of vectors 
in W 1 . Unless otherwise stated, we sum repeated indices over their range, i,j — 
1 , . . . , n where n is the dimension of Ai . We assume that the vector fields and 
1-form densities are defined so that integration by parts gives no contribution at 
the boundary (inclusion of nonzero boundary terms is straightforward). Formula 
( |5~5| ) for ad* c will be useful later. 

By definition, u = (x, v); we will denote u s = x, the spatial part of the phase 
space velocity field. 



The reduced action. We may transform the action (4.1) into the Eulerian de- 
scription as the reduced action 



Srod 



dtl(u, $, $, A, A) 
J dt J dxdvF(x, v, i) Qto|u s | 2 + ^m\u s - v| 2 - g$ + qu s ■ A^j 



1 



+ - / dt j dx 



^ ^ dA 2 



|V x x A| 2 . 

We vary this action with respect to u s , F, $ and A: 



(5.6) 



56 r 



4f (r 



dl I dxilvl F ( 77iu s + m(u s — v) + qA) ■ Su s — q6$ + u s • SA 
SF 

dt I dx[ V x $ 



1 , 
—ml 

2 1 

dA\ 



q& + qu s ■ A 



(5.7) 



dt ) 



Stationary variations in $ and A yield: 



.dA\ 



V x <5$ + 5— J - (V x x A) • (V x x 5 A) 



V x x (V x x A) 



V x • (V x $ + =-qJ dvF(x, v, t), 

q J dvF(x, v,<)u s . 



0_ 

dt 



v x $ 



dA 
~dt 



(5.8) 
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Thus, Maxwell's equations for the electromagnetic field of the plasma are recovered 
by requiring SI — for all variations of the field potentials $ and A. To continue to- 
ward the Euler-Poincare form of the Maxwell- Vlasov equations, one must determine 



the forms of the variations Su s and S F in (5.7). 

According to the general theory, variations in the particle evolution map ip lead 
to variations in the phase space velocity Su of the form 



dw , . dw 
Su = — + [u, wj = — - ad u w . 



(5.9) 



This Euler-Poincare form of the variations may also be verified by a direct tensorial 
calculation, which is given in Holm, Marsden and Ratiu [1997]. The spatial part of 
this equation gives the variation of the spatial part of the field u. 

Variations of the field ip also induce variations of the density F, in the same way 
as the parameter variations are induced in the general theory for the Euler-Poincare 



equations, see equation (5.1). Either from that equation, or by direct calculations, 
these variations are computed to be 



SF = -V. • (Fw), 
which is equivalent to the formula 

S(Fdxdv) = -£ w (Fdxdv). 



(5.10) 



Computation of the variations. With these formulae for Su and SF in place, 
we compute 

5& Ic d — J dt J dudvF ^mu 8 + m(u s - v) + gAj ■ + [u, w] 



- V z • (Fw) ^-(m|u s | 2 + to|u s - v| 2 ) + qu s ■ A - g$ 
Integrating by parts and dropping boundary terms gives 



se 



red — 



dt / (ixrfv W 



d_ 

di 



(Vm(u s + (u s - v) + —A' 

V 771 



ad 



* (Fm(u s + (u s - v) + i-A) 

V 777 



FV y , [ ^m\u s \ 2 + ^rn\u s - v| 2 + qu s ■ A - q$ 



(5.11) 



(5.12) 
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Expanding the ad* term using formula (pM results in 



5& T 



cd 



dt / dxdv w • 



dF ( . . q . 

—m u s + (u s — v) H A 

at V m 



5 , 
Fm-|u s 



(u • V z ) [Fm u s + (u s - v) + —A 

m 



Fm (^u s + (u s - v) 
Fm (u SJ + (u s j - 



-A (V,-u) 

to 



v ? -) + — A 7 - 
m 



■- -/T z ( -m|u s 



We expand the products to obtain 



1 , 



V z u J 
gu s ■ A 



(5.13) 



se 



rod 



J dt J dxdv w | - m + (u s — v) 



<9F 



u • V Z F 



Fm 
f TO 



(u, 



v) + ^A 

m 



u s + (u s - v) 



^A 



q 

-v z $ 



(u sj + (u sj - Vj) + ^-A j J V z //<' + Fmu S jV z u^ + FgA J V z u s j 



+ Fm(u SJ - v J -)V z (u 3 s - vi) + Fqu s:j V z A 3 



(5.14) 



Consider the last two lines of equation (5.14). Upon writing w = (wi,W2), where 
wi,W2 S R 3 , these lines reduce to: 

- Fm(u s + (u s - v)) • (wi • V x + w 2 • V v ) u 
+ Fm(u s + (u s - v)) • (wi • V x + w 2 • V v ) u s 

- FqA • ( Wl • V x + w 2 • V v ) u + FqA ■ (w x • V x + w 2 • V v ) u s 

+ Fqu s • (wi • V x + w 2 • V v ) A - Fm{u s - v) • (w x • V x + w 2 • V v ) v 
= -Fqu s (wi • V x ) A — Fm(u s — v) • w 2 . (5.15) 

The first three lines cancel to zero because they only involve spatial velocity projec- 
tions, where u = u s . The last line follows upon using V x v = and V V A = 0; which 
hold, respectively, because v is an independent coordinate and A is a function of 
space alone. Similarly, and under the additional observation that V z <£ = (V x $, 0) 
because the potential $ also does not depend on velocity, the other three lines of 
equation (5.14) are purely spatial, i.e., the projection onto the last three coordi- 
nates would give zero, and hence the contribution to the variation of the action 
5S re d from w 2 comes only from the calculation in equation (5.15). Stationarity of 
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the action under the velocity components of the variation, W2, then implies: 

Fm(u s — v) = 0, i.e., u s = v. (5.16) 



Consequently, in equation (5.14) we can write u as (v, a) where a is yet to be 
determined, and we can also replace u s — v with zero. On doing this, the contribution 
to the variation of the action from Wi becomes 



S& Te d = I dt I dxdv Wi 



Ov 



OF 

- (mv + qA) ( — + V z • (Fu) 



F m— — h m(v • V x )v + m(a • V v )v 
\ ot 

dA 

q— + gV x $ + qv x (V x x A; 



(5.17) 



Here, we have used standard vector identities in obtaining the result 

w • qF (u sj V z Al - (u • V Z )A) = gFwj • (v x (V x x A)) . (5.18) 



Referring to the continuity equa tion (4.3) for F and using the identities dv/dt = 
and V x v = reduces equation ( 5.17) to: 



66 



red 



Therefore, <56 rc d = implies that 



dt / dxdv wi • F ( ma + gV x $ + q 



dA 



OA 
~dt 



qv x (V x x A) 



ma = - gV x $ - q— + qv x (V x x A). 



(5.19) 



Now consider what the invariance of the Boltzmann function / implies. By 
equation (4.2) and substitution for u = (v, a), we obtain 



-7T7 + V • V x / 

at m 



dA 

V x $ + — ) - v x (V x x A) 



V v / = 0, 



(5.20) 



and so, along with equations ( jO ), we have recovered the full Maxwell- Vlasov system 
from stationarity of the action (5.6) entirely in the Eulerian description. 



6 The Generalized Legendre Transformation. 

Introduction. Before passing to the Hamiltonian description of the Maxwell- 
Vlasov equations, we pause to explain the theoretical background of how one does 
this when there are degeneracies. This section can be skipped if one is willing to 
simply take on faith that one should do the Legendre transformation slowly and 
carefully when there are degeneracies. 

As explained in Marsden and Ratiu [1994], one normally thinks of passing from 
Euler-Poincare equations on a Lie algebra g to Lie-Poisson equations on the dual g* 
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by means of the Legendre transformation. In some situations involving the Euler- 
Poincare equations, one starts with a Lagrangian on q x V* and performs a partial 
Legendre transformation, in the variable £ only, by writing 

SI 

Since 

Sh cm/ 5 *\ / Sl S Z\ C 

^n^rWw = e ' (6 - 2) 



and Sh/8a — —81 /5a, we see that the Euler-Poincare equations (3.4) for £ G 9 and 
a(i) = — a(t)£(t) imply the Hamiltonian semidirect-product Lie-Poisson equations 
for /i e g*. Namely, 

^/j,= -ad*g h/g ^fj,-^oa = {n,h} L p, J^a = -a|^ = {a,h} LP , (6.3) 
with (+) Lie-Poisson bracket on g* x V* given by 

{5 ' /l}iP= -( /i ' ad ™0) + ( a 'l^-Si) ' (6 ' 4) 



If the Legendre transformation (6.1) is invertible, then one can also pass from the 
Lie-Poisson equations to the Euler-Poincare equations together with the equations 
d(t) = -a(t)£(t). 

It is important in this paper to give a detailed explanation that incorporates 
the degeneracy of the parameter dependent system together with the role of sym- 
metry. Unlike the examples considered in Holm, Marsden and Ratiu [1997] such as 
compressible flow or MHD, in the case of the Maxwell- Vlasov system or even the 
Vlasov-Poisson system, the Lagrangian L ao corresponding to the action in equation 
( [5.6] ) is degenerate, since it docs not depend on the variables 3> and v. In other 
words, the degeneracy and corresponding constraints that appear in Vlasov plasmas 
are more serious than for fluids or the heavy top, etc. To deal with this degeneracy, 
we shall use the generalized Legendre transformation in the context of Lagrangian 
submanifolds, as described in Tulczyjew [1977]. This is also related to the Dirac 
theory of constraints (see Dirac [1950]). In particular, we shall take special care to 
ensure that the Hamiltonian formulation of the Maxwell- Vlasov system preserves 
the constraints associated with the degeneracy of its Lagrangian. 



The general construction. Let Q be a manifold and n : T*Q — > Q be the 
cotangent bundle of Q. Then TT*Q is a symplectic manifold with a symplectic 
form that can be written in two distinct ways as the exterior derivative of two 
intrinsic one forms. These two one forms are denoted A and x an d are given in 
coordinates by: 



A = pdq + pdq 



(6.5) 
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and 

X = i>dq~ qdp, (6.6) 

where (q,p) are coordinates for T*Q and (q,P-,q,p) are the corresponding coordi- 
nates for TT*Q. For the intrinsic definitions of these forms, see Tulczyjew [1977]. 

Let L : J — » R be a Lagrangian defined on a submanifold J C TQ called the 
Lagrangian constraint. The Legendre transformation is a procedure to obtain a 
Hamiltonian H : K — * K defined on a submanifold X C T*<2, called the Hamilto- 
nian constraint. The Euler-Lagrange equations are: 

X = dL on J, (6.7) 

while the Hamilton equations are 

X = -dH on K. (6.8) 

The abbreviated expressions ([Tt]) and ( |6.8| ) stand for 

X = d(LoTTr) on (T7r) _1 (J), (6.9) 

and 

X = -d{Ho T ~ 1 ) on (T)- X (if), (6.10) 

where r is the canonical projection r : TT*Q — > T*Q 7 given in coordinates by 
T(q,P,q,p) = (q,p)- The map Tvr is given by TTr(q,p,q,p) = (q,q). 

Both the Euler-Lagrange and Hamilton equations define the same Lagrangian 
submanifold D of TT*Q. The Lagrangian and Hamiltonian L and H are the gen- 
erating functions with respect to the one forms A and \ respectively. 

The generalized Legendre transformation consists of the following steps: 

Step 1. For each {q,p) 6 T*Q define 

d 



K(q,p) = ^(q,q)eT q Q ^ (pq - L(q, q)) = Oj , (6.11) 
and let 

K = {(q,p)eT*Q\K(q,p)^0}. (6.12) 

Assumption. Assume that for each (q,p) G K, the submanifold K(q,p) is con- 
nected. This implies that the stationary value 

st&tq(pq - L(q,q)) (6.13) 

of pq — L(q, q) on K(q,p) is uniquely defined; that is, it does not depend on q. 
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Step 2. Define H : K -> R as follows: 

H(q,p)=st&tq(pq-L(q,q). (6.14) 

The generalized Legendre transformation with parameters and symme- 
try. Now we adapt this methodology to the case of parameter dependent La- 
grangians with symmetry. Let L ao : TG x TQ — > R be a Lagrangian depending on 
a parameter ao <EV* . Assume that G acts on V* on the right and denote by ag the 
action of g G G on a G V*. Assume also the following invariance property: 

Lah(gh, gh, q, q) = L a (g, g, q, q), (6.15) 

for all g, heG, all (g, q) G TQ and all a G V*. A typical element of T*G xT*Q will 
be denoted {g,a g ,q,v q ) or simply (g,a,q,v). For each a G V* and (#, a) G T*G, 
define 

K ao (g,a,q,v) = ^(g,g,q,q) ^r(ag + uq - L a<J (g,g,q,q)^j = 

and ^ ^as- + i/g - L ao (g, g, q, q) ^ = o| . (6.16) 

One can immediately check for any a G V*, h <E G and (g,a,q,v) G T*G x T*Q 
that K aoh (gh, ah, q, v) = K ao (g, a, q, v)h. Define 

K ao = {(g,a,q,u) \ K ao (g,a,q,u) ^ 0}. (6.17) 

Then one can easily prove for any h G G that K ao h — K ao h. Define 

K = {(g, a, q, v, a) \ K a {g, a, q, v) ± 0} . (6.18) 

Then K C T*G x T*Q x V* is an invariant subset under the action of G given by 
(g, a, q, v, a)h = (qh, ah, q, v, ah). Now for each a G V* we define H ao : K ao — > R 

by 

H aa {g,a,q,v) = ag + vq - L aa {g,g,q,q) , (6.19) 

for any (g, g, q, q) G K ao (g, a, q, v). Then, according to the general theory explained 
above, Hamilton's equations are, for each ao G V* , —dH aa = x on K ao , where 

X = adg — gda + vdq — qdv. (6.20) 

One can also easily prove, using the previous equalities, that H ao (g, a, q, v) has the 
following invariance property, 

H aoh (gh,ah,q,v) = H ao (g,a,q,v). (6.21) 

Let s* be the dual of the semidirect product Lie algebra s = Q®V . Then define 
K C s* x T*Q by 

JC = {(a, q, v, a) G S* x T* Q \ (e, a, q, v, a) G K} , 



LAUR-97-3326 



Maxwell- Vlasov system in Euler-Poincare form 



19 



and the Hamiltonian hjc ■ IC — > K by hjc{a,a,q,v) = H a {e,a,q,v). Thus, h/c is 
the restriction to IC C s* of the right invariant Hamiltonian H : K — > K given 
by -£f(<?, a, 9, f, a) = H a (g, a, q, v). Then, by a natural generalization of semidirect 
product theory to include constrained Hamiltonian systems, we have that Hamil- 
ton's equations on JC C s* generated by hjc give the evolution of the system on 
IC determined by the Poisson-Hamilton equations / = {/, hjc} on the Poisson sub- 
manifold IC C s* x T* Q, where the Poisson structure is defined in a natural way. 
More precisely, we have the Dirac brackets on K (see for instance Dirac [1950] or 
Marsden and Ratiu [1994]) which, by reduction, give the brackets on IC. This is 
the abstract procedure underlying the computations we do in the specific case of 
plasmas given in the next section. 



7 Hamiltonian formulation 

We now pas s to the corresponding Hamiltonian formulation of the Maxwell- Vlasov 
system (|2.4|)-(|2.lij) in the Eulerian description by taking the Legendre transform of 



the reduced action (5.6) 



The role of the general theory. From the geometrical point of view, we simply 
apply the generalized Legendre transformation described abstractly in §^ to the 
degenerate Lagrangian 

^/o^oO.^*.*. A, A) :TT xTV -xTA^M.. 

This Lagrangian is degenerate because it does not depend on the variables 4> and 
v. The theory described in ^ may be applied to this action on T(T x V X A). 
The action of the group Diff (TM 3 ) on the factor T for this Lagrangian is given as 
before, while the actions on the factors V and A are trivial. It is easy to see that 
the Hamiltonian constraint for each / is Kf C T*(Diff(TR 3 ) x V x A), defined by 
the conditions 

5L , 5L 

* = — = and m„ = — = 0. 

These conditions impose constraints, which for consistency must be dynamically 
preserved. 

Calculation of the transformed equations. We will perform the calcula- 
tions in detail, working with the reduced Lagrangian rather than the Lagrangian 
L/ 0iVo (^,V>,$,$, A, A) and setting ip n = 1 as usual. 

We start with the action ( |5.6| ) for the Maxwell- Vlasov system in the Eulerian 
description, 

6 rcd (u,$,$,A, A) = 

Jdtjdxd-v F(x,v,<) Qm|u s | 2 + ^m|u s - v| 2 - 0$ + qu s ■ A 



dt I dx 



dA 



2 



|V x xA| 2 . (7.1) 
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This leads immediately to 

61 „ x dA „ ffrn . 
-r = V x $ + — = -E, 7.2 
SA dt 

and so (minus) the electric field variable E is the field momentum density canonically 
conjugate to the magnetic potential. Let us define the material momentum density 
in six dimensions, 

m-|. (7.3) 

We write m = (m s ,m„), where m s is the projection of m onto the first three 
cordinate positions, and m„ is the projection onto the last three places. We think 
of m s and m„ also as vectors in six dimensions. From the Lagrangian we see that 



m s = Fl mu s + m(u s — v) + qA\ and m„ = 0. (7-4) 



Proceeding with the Legendre transform of our action (7.1) results in a corre- 
sponding (reduced) Hamiltonian function written in terms of the velocities, 

dxdv F(m|u s | 2 - -rn|v| 2 + g<f> 



+ i J dx(|E| 2 + |V X x A| 2 + 2E • V x $) , (7.5) 

where a denotes the projection of u onto its last three entries. Transforming this 
to the momentum variables gives 

f 1 1 

h = I dxdv |m s + mFv - qFA\ 2 mF\v\ 2 + qF§ + m„ • a 

J AFm 2 

+ i J dx(|E| 2 + |V X x A| 2 + 2E • V x $) . (7.6) 
The variation of this Hamiltonian with respect to m, a, E, A, F and $ is given by 



Sh = dxdv 



u • 6m + m v ■ 6a — qFu s ■ 6 A + qF5Q 
im|u s | 2 + im|u s - v| 2 + qu s ■ A - q$J 6F 

+ J dx(E + V x $) • 6E - (V x • E)5$ + V x x (V x x A) • SA. (7.7) 

This expression allows one to read off the evolution equations for the electromagnetic 
field: 

dA Sh dA 

lH=-SV = - E - V ^ Le - E = - Vx$ -^< 

Sh f f 

— = = -V x ■ E + q / dvF, i.e., V x ■ E = q / dvF := p, 

— = — = V x x (V x x A) - q / dvFu Sl i.e., — = V x x B - j . (7.8) 
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Note that the constraint 5h/5Q = (Gauss' law) arises from the absence of $ 
dependence in I. 

The general theory of §^ shows that F is an element of the second factor of the 
semidirect product and so its evolution is given by Lie dragging as a density. Like- 
wise, / is Lie dragged as a scalar and rrii satisfies a Lie-Poisson evolution equation: 

dF 

— = -V.-(Fu), 



dt 
dj_ 
dt 

dm; d 4 d „• „ d Sh 



V z /, 



.miU 1 — m;——v? — F (7.9) 

dt dzi 1 dz l dz l 5F v ; 

The first two of these equations reflect the assumptions that were made in the 
definitions of / and F, whilst the last equation encodes the dynamics for the system. 
We first consider the case where the momentum component i takes the values 4, 5, 6. 
In this case, 

® i 9 ,,- _ d (\ . l2 1 l2 . , 
= ^sj-^-rU J +m vjT —u ■> - F T —(-m\u s \ + -ra|u s -v| +qu s ■ A - q$ 



dt J dz l J dz l dz l \2 

3 Qu^ 

dkj | p d 
dz l dz 



1 F *s 3 ^d+F q — $, (7.10) 



where i = 4, 5, 6. In the second line of equation ( 7.10| ), we have substituted for m. 



from equation (7.4) and rearranged terms. Here m„ = 0, because I does not depend 
on v. Setting m„ = initially in equation ( 7.10| ) ensures that m v = persists 



throughout the ensuing motion; for potentials <I> and A that are independent of v, 
and provided the constraint holds that u s — v, as in equation Q5.16 ). Likewise, 



the Gauss' law constraint imposed by 5h/5<S> = also persists during the ensuing 
motion, as seen from the last equation of ( |7.8| ) and the first equation of (|7.9|), 
provided the constraint u s = v holds and F vanishes in the limit as |v| — > oo . 

The spatial part of the evolution equation of m will produce the required single- 
particle dynamics. From equation (|7.9|), we have 



dim d d d Sh 

~~dT dzl miU ~ m ^ u ~ F d?SF- (7 ' n) 



Setting i — 1,2,3, in equation ( 7.11 ), then substituting for Sh/SF and using the 
relations 

m s = F(mu s + m(u s - v) + qA), m„ = 0, 

and 



V v $ = 0, V V A = 0, 
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yields the spatial components of the motion equation, 

dm si d d ■ 

dt dzi J dz l 

- F^- Qm|u s | 2 + im|u s - v| 2 + qu s ■ A - q$j . (7.12) 

Substituting for m s and then using the continuity relation dF/ dt + V z • (Fu) = 
gives 

du si dA st d , d did* 

m tt: \~ Q 7rr~ = - u J -—mu si - u J s -—qA sl - qAj — -u J - q— 

dt dt dz J dzi oz l az 

d , .,19 



2 



+ ? 9?( U8 ' A )-2 m 9? |Us_v| - (7 ' 13) 
Rearranging this equation results in 

{ d Q \ 1 d 

m \ + uj lTl = 9 E + <?u s x (V x x A) - -to-— |u s - v| 2 . (7.14) 



dt dzi) s H H * v x ' 2 dz v 

We may now evaluate this on the constraint set u s = v and thereby obtain the 
Lorentz force, 

ma = q(E + v x B), (7.15) 

where a is the acceleration of a fluid parcel (the last three components in u = va). 
As we have seen, in this Hamiltonian formulation of the Maxwell- Vlasov equations 
in the Eulerian description, the acceleration a in u is a vector Lagrange multiplier 



which imposes m v = 0. Equation (7.15) provides an expression for this Lagrange 
multiplier in terms of known dynamical variables and as a consequence we regain 
the equation for the acceleration of a charged particle in an electromagnetic field. 
The momentum constraint m,, = remains invariant when the electromagnetic 
potentials are independent of the phase space velocity coordinate v and the velocity 
constraint u s = v holds. Perhaps not unexpectedly, one finds that V z • u = 0. 
Also, (minus) the electric field is canonically conjugate to the vector potential, and 
the electrostatic potential $ plays the role of a Lagrange multiplier which imposes 
Gauss's law. Thus, our Hamiltonian formulation augments the usual Maxwell- 
Vlasov description of plasma dynamics by self-consistently deriving the particle 
acceleration by the Lorentz force ma. = q(E + v x B) instead of assuming it a 
priori. 



The Poisson Hamiltonian structure. The general theory outlined briefly in 
also leads to the Poisson bracket structure for the Maxwell- Vlasov theory on 
the Hamiltonian side. However, our Hamiltonian description has a redundancy, 
namely the information for the particle trajectories can be recovered from the spatial 
plasma density. Explicitly, if we let H(f) = (l/2)|v| 2 + $/(x) be the single particle 
Hamiltonian determined by the plasma density /, then the flow of this Hamiltonian 
function can be identified with the particle evolution map ip- We can also think of 
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this as a constraint on the level of equations of motion, as the Hamiltonian vector 
field of H(f) must equal the time derivative of the map ip, i.e., the particle velocity 
field in phase space. In other words, as is well known, the particle dynamics is 
completely determined by the plasma density dynamics. This may be regarded as 
a constraint on the system that leads to the elimination of the forward map as a 
dynamical variable. This "redundancy" is of course one of the sources of degeneracy 
of the Lagrangian and Hamiltonian structures. 

Thus, the constraint explicitly enforcing this consistency condition leads to a 
further "reduction" which again may be handled by the Dirac theory of constraints 
to arrive at the Hamiltonian structure in terms of the variables F (or equivalently / 
in view of the canonical nature of the particle transformations) and the electromag- 
netic potentials. The resulting Poisson bracket structure is given by the Lie-Poisson 
structure for the /'s plus the canonical structure for the electromagnetic potentials, 
which was the starting point for Marsden and Weinstein [1982], who carried out the 
reduction of this bracket with respect to the action of the electromagnetic gauge 
group to obtain the final Maxwell- Vlasov bracket on the space with variables /, E 
and B. This procedure was motivated by and corrected a bracket found by ad hoc 
methods in Morrison [1980]. We need not repeat this construction. 



The Kelvin-Noether theorem. A final result worth mentioning is a Kelvin's 
theorem for the Maxwell- Vlasov particle dynamics. These dynamics, given in the 
last equation in (7.9), may be rewritten as 

^ + ^„)(im i (fo i )+d^=0, (7.16) 
so that 

4/ i^^O, (7.17) 

dtJ l{t) F V ' 

for a loop j(t) which follows the particle trajectories in phase space. The Kelvin 
circulation integral in phase space, 

/=/ ^ mi dz\ (7.18) 

J-y(t) * 

may be evaluated on the invariant constraint manifold = as 



I = f (rmisi + qAi)dx l . (7.19) 
i 7 (t) 

We recognize this integral as the Poincare invariant for the single particle motion 
in phase space. 

The above result follows from the abstract Kelvin-Noether theorem by letting 
C := {7 : S 1 -> TR 3 | 7 continuous} be the space of continuous loops in single 
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particle velocity phase space and letting the group Diff(TR 3 ) act on C on the right 
by (Vi 7) S Diff(T]R 3 ) x C i-» 7 077 eC. The quantity K, is chosen to be 

{K,{ 1 ,F),a) = jja. (7.20) 

The abstract Kelvin-Noether theorem for the Maxwell- Vlasov equations in Euler- 
Poincare form then reproduces the version of Kelvin's theorem given in (7.17). 



8 Conclusion 

In this paper we have cast Low's mixed Eulerian-Lagrangian action principle for 
Maxwell- Vlasov theory into a purely Eulerian description. In this description we 
find that Maxwell- Vlasov dynamics are governed by the Euler-Poincare equations 
for right invariant motion on the diffeomorphism group of K™ (n = 6 for three di- 
mensional Maxwell- Vlasov motion). These equations were recently discovered by 
Holm, Marsden and Ratiu [1997] who investigated the class of Hamilton's principles 
which are right invariant under the subgroup of the diffeomorphisms which leaves 
invariant a set T of tensor fields in the Eulerian variables. The Maxwell- Vlasov 
motions invariant under this subgroup are the steady Eulerian solutions, which, 
thus, arc identified as relative equilibria. This identification of steady Eulerian 
Maxwell- Vlasov solutions as right invariant equilibria places these solutions into 
the Hamiltonian framework required for investigating their nonlinear stability char- 
acteristics using, e.g., the energy-Casimir method (see Holm et al. [1985]). It was 
this stated goal that first motivated Low to write his Lagrangian for Maxwell- Vlasov 
dynamics. 

Thus, our formulation of a purely Eulerian action principle and its associated 
Euler-Poincare equations and Hamiltonian framework advances Low's original in- 
tention of using his action principle for studying stability of plasma equilibria by 
placing the entire Maxwell- Vlasov equations (including the particle dynamics, field 
dynamics and probability distribution dynamics) into one self-consistent Hamilto- 
nian picture in the Eulerian description. (As we discussed, Low used mixed aspects 
of both Eulerian and Lagrangian phase space descriptions in his action principle.) 

Our Eulerian Hamilton's principle for Maxwell- Vlasov dynamics is constrained, 
and all of the corresponding Lagrange multipliers have been resolved. This Hamil- 
ton's principle is thus available for further approximations, e.g., by Hamilton's prin- 
ciple asymptotics (see, e.g., Holm [1996]). 

In summary, we have taken an existing action, due to Low [1958], for the 
Maxwell- Vlasov system of equations and demonstrated how to rederive this system 
as Euler-Poincare equations. The Euler-Poincare form emerges from Hamilton's 
principle for a system whose configuration space is a group and whose action is 
right invariant under a subgroup. This situation commonly appears in the Eule- 
rian description of continuum mechanics. In the case of continuum mechanics, the 
dynamics takes place on the group of diffeomorphisms and the Eulerian variables 
are invariant under a subgroup of the diffeomorphism group. (This subgroup corre- 
sponds to steady Eulerian flows with non-zero velocity and vorticity.) We showed 
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that this situation also occurs for the Maxwell- Vlasov equations of plasma dynam- 
ics in the Eulerian description, by showing that the variations considered take the 
appropriate form, and then deriving the Maxwell- Vlasov equations from the Hamil- 
ton's principle for the right invariant action (5.6) in Eulerian variables. We then 
passed to the Hamiltonian formulation of this system and found its Lie-Poisson 
structure. 

As discussed in the introduction, the Euler-Poincare form of the dynamics is 
naturally adapted for applying Lagrange-D'Alembert methods for geometrical con- 
straints and control as in Bloch et al. [1996]. In future work, our Euler-Poincare 
form of the Maxwell- Vlasov system shall be implemented to describe the control 
features of a plasma driven by an external antenna, following the lines of inquiry 
begun in the oscillation center approximation for plasmas by Similon et al. [1986]. 
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